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ABSTRACT 

The monopole confinement mechanism in the abelian projection 
of lattice gluodynamics is reviewed. The main topics are: the abelian 
projection on the lattice and in the continuum, a numerical study 
of the abelian monopoles in the lattice gauge theory. Additionally, 
we briefly review the notation of differential forms, duality, and the 
BKT transformation in the lattice gauge theories. 



1 Lectures given by Mikhail Polikarpov at the Workshop "Confinement, Duality and Non-Perturbative 
Aspects of QCD", Cambridge (UK), 24 June - 4 July 1997. 



1 Introduction 



In these lectures we give an introduction to the theory of the confinement of color in lattice 
gauge theories. For the sake of self-consistency, we explain all definitions. The reader is 
supposed to be familiar with basic notions of lattice gauge theory. 

Lattice field theories were originally formulated U in order to explain the confinement of 
color in nonabelian gauge theories. The leading term of the strong coupling (l/g 2 ) expansion 
in lattice gauge theories yields the area law for the Wilson loop. The D dimensional theory 
is reduced to the 2-dimensional one, the field strength tensors are independent on different 
plaquettes and the confinement has a stochastic nature. The expectation value of the Wilson 
loop is simply < P > s , where < P > is the expectation value of the plaquette, S is the area 
of the minimal surface spanned on the loop. Therefore, the string tension is a = — In < P >. 

The numerical study of lattice gauge theories initiated by Creutz || shows that the strong 
coupling expansion of lattice gauge theory has no relation to continuum physics: the effects 
of lattice regularization are very strong, the results are not rotationally invariant, and there 
is no scaling for hadron masses^]. 

In the weak coupling region (where the existence of the continuum limit was found in 
numerical calculations), the confinement mechanism has no stochastic origin. There are 
several approaches to explain color confinement. We will describe the most traditional one, 
namely, confinement caused by the dual Meissner effect. The linear confining potential can 
be explained by the formation of a string (flux tube) connecting a quark and an anti-quark. 
A well-known example of the string-like solution of the classical equations of motion is given 
in Ref. []3J. If we have a medium of condensed electric charges (a superconductor), then 
between the monopole and anti-monopole an Abrikosov string is formed, see Figure |l|(a). To 
explain the confinement of electric charges, we need a condensate of magnetic monopoles, 
Figure |](b). This simple qualitative idea was suggested by 't Hooft and Mandelstam ||]. 
There is a long way from this picture to real QCD. First, we have to explain how the abelian 
gauge field with monopoles is obtained from the non-abelian gauge field. Secondly, we have 
to explain why a dual superconductor is involved here. We discuss these two questions in 
these lectures. 

At present, we have no analytic proof of the existence of the condensate of abelian mag- 
netic monopoles in gluodynamics and in chromodynamics. However, in all theories allowing 
for an analytical proof of confinement, the latter is due to the condensation of monopoles. 
These analytical examples are: compact electrodynamics M, the Georgi-Glashow model 0, 
and super-symmetric Yang-Mills theory 0. 

On the other hand, many numerical facts (some of these are discussed in these lectures) 
suggest that the vacuum in SU{2) and SU(3) lattice gauge theories behaves as a dual 
superconductor. As an illustration, we give two figures obtained by numerical calculations 
in SU (2) lattice gluodynamics. 

In Figure 0, taken from Ref. ||, the action density (vertical axis) of the SU(2) fields is 
shown. The two peaks correspond to the quark-anti-quark pair, the formation of the flux 

"Here, scaling means independence of the ratios of the hadron masses mk/rrii on the parameters of the 
theory such as the bare coupling and the cutoff parameter. 
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(a) (b) 

Figure 1: The Abrikosov string between the monopoles in the superconductor (a) and an 
analogue of the Abrikosov string between the electrically charged particles in the dual su- 
perconductor (b). 




Figure 2: The action density distribution || at /3 = 2.635, the distance between quark and 
anti-quark sources is ~ 1.35/m. Two horizontal axes correspond to two spatial lattice axes. 

tube is clearly seen. In Figure |3], taken from Ref. M, the abelian monopole currents near 
the center of the flux tube formed by the quark-anti-quark pair are shown. It is seen that 
the monopoles wind around the center of the flux tube just as the Cooper pairs wind around 
the center of the Abrikosov string. 

We first explain how to get the abelian fields and monopoles from the non-abelian fields. 
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Then we present the results of numerical studies of the confinement mechanism in lattice 
gluodynamics. 

All technical details such as the formalism of differential forms on the lattice are given 
in the Appendices. 

2 Abelian Monopoles from Non-Abelian Gauge Fields 

In this Section we discuss the question how to obtain the abelian monopoles from non-abelian 
gauge fields. 

2.1 The Method of Abelian Projection 

The abelian monopoles arise from non-abelian gauge fields as a result of the abelian pro- 
jection suggested by 't Hooft [|K|. The abelian projection is a partial gauge fixing under 
which the abelian degrees of freedom remain unfixed. For example, the abelian projection 
of a theory with SU(N) gauge symmetry leads to a theory with [[/(l)]^ -1 gauge symmetry. 

Since the original SU(N) gauge symmetry group is compact, the remaining abelian gauge 
group is also compact. But the abelian gauge theories with compact gauge symmetry group 
possess abelian monopoles. Therefore SU (N) gauge theory in the abelian gauge has abelian 
monopoles. 
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Figure 3: The monopole currents around the string tube which is formed between the static 
quark and anti-quark in gluodynamics, Ref. ||. 
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First, consider the simplest example of F± 2 abelian projection for SU(2) gauge theory. 
This gauge is defined by the following condition: 



F\ 2 (x) = diagonal matrix . (1) 

It is easy to fix to this gauge, since under gauge transformations the field strength tensor 
F12 transforms as 

F 12 (x) -> F[ 2 (x) = Q + (x)F 12 (x)Q(x) . (2) 

If we fix to the F 12 gauge, then the field strength tensor F[ 2 (x) is invariant under U(l) gauge 
transformations: 

F[ 2 {x) = % {1) {x)F[ 2 {x)n u{1) {x) , (3) 

where 

Vu(i){x) ={ e e -<a(*) J > a E t°' 27r ) • ( 4 ) 

Therefore, the gauge condition (0) fixes the SU(2) gauge group up to the diagonal U(l) 
subgroup. 

The SU (2) gauge field 

Ml -V <«> 

transforms under the gauge transformations as 

-> Q+A^Q - -Q+d^Q . (6) 

If we fix to the F\ 2 gauge, then under the remaining U(l) gauge transformation (Q) the 
components of the nonabelian gauge field A transform as 
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<u,-<u,--d lt a, o M = 4; A+ ^ e 2ia A+ , Al = A\ + iA\. (7) 
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Thus, in the abelian gauge the field plays the role of the abelian gauge field and the field 
A^ L is a charge 2 abelian vector matter field. 

We thus obtain abelian fields from non-abelian ones. It occurs that abelian projection 
is also responsible for the appearance of abelian monopoles. If a is a regular abelian field, 
then divif = 0, since H = curia. But the abelian gauge field may be non-regular, since 
the matrix of the gauge transformation Q may contain singularities. The nonabelian field 
strength tensor F^ v is not invariant under singular gauge transformations: 

F» U [A] - F,M {n) ] = V+F^AP + FSfM , (8) 
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where 

F s ;m = --n + (x)[d^ u - a v d^nx) . (9) 

Therefore, if we fix to the abelian gauge (|TJ), using the singular gauge rotation matrices, 
the abelian field strength tensor 

Uu = d^a,, - d„a^ (10) 
may contain singularities (the Dirac strings) 

fixv = fju, + f% , Sfj.ua/3dufal3 — anci £liva(3d v fa(3 7^ , (11) 

and therefore, div/7 7^ 0. 

The charge of the monopole can be calculated [1]J by means of the Gauss law. Let us 
choose an abelian monopole in a certain time slice and surround it by an infinitesimally 
small sphere S. The monopole charge is 

m = -fHdB=—J do^ e^ va(3 f a(3 = J do^ e^ va(3 Tr [Q + d a QQ + dpQ] 

= --^ fda^e^d^rin+d^} = 0, ±-, . . . . (12) 

The quantization of the abelian monopole charge is due to topological reasons: the surface 
integral in eq. (|12D is equal to the winding number of SU(2) over the sphere S surrounding 
the monopole. The physics of quantization is simple: the electric charge is fixed by the gauge 
transformation (|?p, and magnetic charge should obey the Dirac quantization condition. The 
last integral can be seen as the magnetic flux through the Dirac string that arises as a 
consequence of the gauge singularity. 



2.2 Various Abelian Projections 

There is an infinite number of abelian projections. In the previous section we have considered 
the F12 abelian gauge. Instead of the diagonalization of the tensor component F12 by the 
gauge transformation, we can diagonalize any operator X which transforms under the gauge 
rotation as follows: X — > Q + XQ. Each operator X defines an abelian projection. At finite 
temperature one can consider the so-called Polyakov abelian gauge which is defined by the 
diagonalization of the Polyakov line. 

The most interesting numerical results are those obtained in the Maximal Abelian (MaA) 
gauge. This gauge is defined by the maximization of the functional 

max R[A^] , R[A] = - J d 4 x [(^) 2 + (Aj) 2 ] , (13) 

The condition of a local extremum of the functional R is 

(^±^4)^ = 0. (14) 

Clearly, this condition (as well as the functional R[A}) is invariant under the U(l) gauge 
transformations (0). The meaning of the MaA gauge is simple: by gauge transformations 
we make the field A M as diagonal as possible. 
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2.3 Abelian Projection on the Lattice 



The SU(2) gauge fields JJ\ on the lattice are defined by SU{2) matrices attached to the 
links I. These lattice fields are related to the continuum SU(2) fields A: U X:fl = e la9A ^ x \ 
here a is the lattice spacing. Under the gauge transformation, the field Ui transforms as 
Ux,n = Ux,^ &x+fa the matrices of the gauge transformation are attached to the sites x of 
the lattice. 



The Maximal Abelian gauge is defined on the lattice by the following condition 11 



m^R[U?), R[U l }=Y J Tr[ozU+ozU l ], 1 = {x, ^ . (15) 

i 

This gauge condition corresponds to an abelian gauge, since R is invariant under the gauge 
transformations defined by the matrices (§]). 

Let us parametrize the link matrix U in the standard way 

ui={ coscpie Z xl sin ^;), (i6) 

V — sm ipie Xl cos fie 1 ) 

where 9, \ £ [ — + 71 ") and (p G [0, it). In this parameterization, 

R[Ut] = $>os2<^. (17) 
i 

Thus, the maximization of R corresponds to the maximization of the diagonal elements of 
the link matrix (|l6l). 

Under the U(l) gauge transformations, the components of the gauge field (|16D are trans- 
formed as 



®x,[i * ®x,[i Q-x (%x+jl j Xx,/i * Xx,/i ®-x (%x+jl i ^Px,p, * (1^) 

Therefore, in the MaA the gauge, the field 9 is the U (1) gauge field, the field x is the abelian 
charge 2 vector matter field, the field <p is the non-charged vector matter field. 



2.4 Monopoles on the Lattice 

A configuration of abelian gauge fields 9i can contain monopoles. The position of the 
monopoles is defined by the lattice analogue of the Gauss theorem. Consider the elementary 
three-dimensional cube C (Figure |(a)) on the lattice. 

The abelian magnetic flux H through the surface of the cube C is given by the formula 

Z7r P&dC 

where 9p is the magnetic field defined as follows. Consider the plaquette angle 9p = 9\+9 2 — 
9s — #4 = d9, the 9iS are attached to the links % which form the boundary of the plaquette 
P, Figure f§(b). The definition of 9p is 9p = 9p + 2nk, where the integer k is such that 
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(a) (b) 

Figure 4: The magnetic flux through the boundary of the cube C (a) and lattice field strength 
tensor (b). 

—7i < 9p < 7i. The restriction of 9p to the interval (— 7r,7r] is natural since (as we point 
out in Appendix A) the abelian action for the compact fields 9\ is a periodic function oiOp. 
Equation ([19]) is the lattice analogue of the continuum formula m = § HdS. Due to the 
compactness of the lattice field 9 (— n < 9 < tt) there exist singularities (Dirac strings), and 
therefore, divif ^ 0. 




Figure 5: The definition of the magnetic monopole current. 

The magnetic charge m defined by eq. fll~9|) has the following properties: 

1. m is quantized: m = 0, ±1, ±2; 

2. If m 0, then there exists a magnetic current j. This current is attached to the link 
dual to cube C, Figure [5|. 



3. Monopole currents j are conserved: 8j = 0, the currents form closed loops on the AD 
lattice. The proof is given in Example 3 of Appendix A. 



2.5 Vacuum of Gluodynamics As a (Dual) Superconductor 

As we have just shown, the non-abelian gauge field in the abelian projection is reduced 
to abelian fields and abelian monopoles. The static quark-anti-quark pair in the abelian 
projection becomes the electric charge-ant i- charge pair, and if the monopoles are condensed, 
then the quark and anti-quark are connected to each other by an analogue of the Abrikosov 
string. A more detailed discussion of the confinement in the abelian projection is given 



in Refs. [ID . Thus, in order to justify the monopole confinement mechanism we have to 
prove the existence of the monopole condensate. For lattice gluodynamics we have a lot 
of numerical facts which confirm the monopole confinement mechanism. We discuss these 
at the end of this section. But first we give an analytical example which shows how the 
monopole condensate appears in compact electrodynamics. 

2.5.1 Compact £7(1) Lattice Gauge Theory 

We show that lattice compact electrodynamics can be represented as a dual Abelian Higgs 
model, with the Higgs particles being monopoles. The partition function of the compact 
£7(1) gauge theory can be written as 

Z = J V6i exp{-S{d6)} , (20) 

— 7T 

+ 7T +7T 

where, as below, J T>6[ = 11/ d$z is the integral over all link variables. In the continuum 

limit (o-O), S(ep)<xf* v y* 

By means of the duality transformation (see Appendix B), the theory (|20| ) can be rewrit- 
ten in the form 

Z = const. ■ exp{-S d (d*n)} , (21) 

*n(*C 1 )G^ 
<5*n=0 

where d*n = *n\ + *ri2 — *n 3 — *n 4 is the plaquette constructed from the integers n (see 
Figure |6]). The dual action is defined by eq. (|B.3| ). 

The integer valued variable n X)fl is dual to the field 9 XtlM . In the continuum limit, 8 X)fl oc 
A^(x) is the ordinary £7(1) gauge field and n x ^ oc B^ix) is the dual £7(1) gauge fieldQ, 

9fiA u — 8 U A^ = T^E^uapdaBp. 

The partition function (|2l| ) can be represented as that of the (dual) Abelian Higgs model. 



+oo 

Z = const. ■ Hm Km / V*B J exp{-S AH (*B, *$)} , (22) 



"Tii the continuum notations we do not use the symbol "*" to denote the dual fields. 
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n 2 





Figure 6: Field strength tensor on the dual lattice, 
where the action of the (dual) Abelian Higgs model is 

S AH (*B, *$) = £ ^(d*5) + I ^ ^ - e^-*^/ + A E(l*^-| 2 " 7) 2 • (23) 

P * X fi=l 

Here d*B is the plaquette variable constructed from the link variable *Bi, where *Bi is the 
dual gauge field. In the continuum limit, we have S(d*B) — > f3(d fJ ,B v — d u B^) 2 . The second 
term in eq.(|23|) has the following continuum limit: \$ x — e l * Bx ^*$> x +p\ — > {(d^ — iB^)§\ 2 . 
The action S AH is invariant under the following gauge transformations: *B — > *S — d*a, 
*$ _> *$e^ a . 

In the London limit, A — > oo, the mass of the Higgs particle tends to infinity. The radius 
of the Higgs field |*$| is fixed and only the phase *cp of the Higgs field *$ = l^le 4 *^ is a 
physical degree of freedom. In the unitary gauge *(p = 0, the action of the Abelian Higgs 
model (@) is 

S AH = E S d {d*B) + - E E I 1 ~ e l * B ^f . (24) 
p 2 x ^ =4 

In the limit 7 — > 00, the photon mass = 7//? becomes infinite and the field is 
equal to zero modulo 2n: *B X ^ = 2ii*n x ^. In this limit, the partition function ( p4|) reduces 
to 

Z = const. ■ E e- sd{d * n) , S(d*n) = 4vr 2 /3 (d*n) 2 . (25) 

5*n=0 

Therefore, the compact U(l) gauge theory is equivalent to the dual abelian Higgs model 
in the double limit 

Compact U(l) gauge theory = lim lim (Abelian Higgs model) (26) 

7^00 A— >oo 

The gauge fields B\ in eq. (|23|) are dual to the original gauge fields 9i, and these inter- 
act via the covariant derivative with the Higgs field $. Therefore, the field $ carries the 
magnetic charge, and due to the Higgs potential in eq.(^3|), these monopoles are condensed 
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at the classical level. It is well known that in the quantum AD compact electrodynamics 
there exists a confinement-deconfinement phase transition. It can be shown by numerical 



calculations [Eg, [13[ that in the confinement phase the monopoles are condensed and that in 



the deconfinement phase the monopoles are not condensed. 
2.5.2 What is the Theory Dual to Gluodynamics? 

In any abelian projection, lattice gluodynamics corresponds to some abelian gauge theory. 
This abelian theory, in general, is very complicated and non-local. Nevertheless, in the 
Maximal Abelian projection, numerical calculations show that the abelian monopoles are 
important degrees of freedom, and that they are responsible for the confinement. Moreover, 
as we show in the next section (see also Refs. |T^, the distribution of monopole currents 
indicates that, at large distances, gluodynamics is equivalent to the dual Abelian Higgs 
model, the Higgs particles are abelian monopoles and these are condensed in the confinement 
phase. 

3 Numerical Facts in AD SU(2) Gluodynamics 

The standard scheme of numerical calculations can be described as follows. 

1. Generate the lattice Yang-Mills fields Ui using the standard Monte-Carlo method. 
Thus we obtain the configurations of fields distributed according to the Boltzmann 
factor e~ s 



-•YM 



2. Perform the abelian gauge fixing and extract abelian gauge fields from the non- abelian 
ones. In the case of the MaA projection, the abelian gauge fixing is a rather time- 
consuming problem. 

3. Extract abelian monopole currents from the abelian gauge fields. As mentioned above, 
the monopole currents form closed paths on the dual lattice. In Figure [7] we show the 
abelian monopole currents for the confinement (a) and the deconfinement (b) phases. 
It is seen that in the confinement phase the monopoles form a dense cluster, and there is 
a number of small mutually disjoint clusters. In the deconfinement phase the monopole 
currents are dilute. 

4. Calculate expectation values of various operators using the monopole currents. Below 
we discuss a number of numerical facts which show that abelian monopoles in the MaA 
gauge are the appropriate degrees of freedom to describe confinement. 



3.1 Fact 1: Abelian and Monopole Dominance 

The notion of the "abelian dominance" introduced in Ref. means that the expectation 
value of a physical quantity < X > in the nonabelian theory coincides with (or is very 
close to) the expectation value of the corresponding abelian operator in the abelian theory 
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obtained by abelian projection. Monopole dominance means that the same quantity can be 
calculated in terms of the monopole currents extracted from the abelian fields. If we have 
N configurations of nonabelian fields on the lattice, the abelian and monopole dominance 
means that 

1-1 1 

^ £ XiUnanaOelian) = ^E X'(U abelian ) = - £ X"(j) . (27) 
con} con} con} 

Here each sum is taken over all configurations; Uabeiian — is the abelian part of the non- 
abelian field U nonabelian-, j is the monopole current extracted from U a ij e u an . It is clear that 
jr Hconf X(V r n onabeiian) is a gauge invariant quantity, while the abelian and the monopole con- 
tributions depend on the type of abelian projection. In numerical calculations the equalities 
(P7|) can be satisfied only approximately. 

Among the well-studied problems is that of the abelian and the monopole dominance for 
the string tension fL6|, |17|, [Ej| p]. In this case, X {U nonabelian) = <rsu(2), ^{U abelian) = (Tuq.) 



and the string tension asu(2) i a u{i)) is calculated by means of the nonabelian (abelian) 
Wilson loops, TrYl ieC Ui ( Iliec* 6 ^)- An accurate numerical study of the MA projection of 
SU(2) gluodynamics on the 32 4 lattice at (5 — 2.5115 is performed in Ref. JTP| . 



The corresponding string tension can be taken from the potential between a heavy quark 
and antiquark: V(r) = Vo + crr — -, where r is the distance between the quark and antiquark. 



12 



DC 
> 




Figure 8: Abelian and nonabelian potentials (with the self energy Vq subtracted), Ref. [19 
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Figure 9: The abelian potential (diamonds) in comparison with the photon contribution 
(squares), the monopole contribution (crosses) and the sum of these two parts (triangles), 
Ref. P||. 



The abelian and the nonabelian potentials are shown in Figure |[ The contribution of the 
photon and the monopole parts to the abelian potential is shown in Figure 

The differences in the slopes of the linear part of the potentials in Figure |8| and Figure ^| 
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yield the following relations: &u(i) ~ 92% <Jsu(2), &j ~ 95%au(i), where Oj is the monopole 
current contribution to the string tension. It is important to study a widely discussed idea 
that in the continuum limit (f3 — > oo), the abelian and the monopole dominance is exact 

a SU(2) = = CTj. 

There are many examples of abelian and monopole dominance in the MA projection. The 
monopole dominance for the string tension has been found for the SU(2) positive plaquette 
model in which Z 2 monopoles are suppressed |20 | , and also for the SU(2) string tension at 



finite temperature [ETJ], and for the string tension in SU(3) gluodynamics Abelian and 
monopole dominance for SU(2) gluodynamics has been found in |]23| , p4| for the Polyakov 
line and for critical exponents for the Polyakov line, for the value of the quark condensate, 
for the topological susceptibility and also for the hadron masses in quenched SU(3) QCD 
with Wilson fermions [E5fl. 



3.2 Fact 2: London Equation for Monopole Currents 

In the ordinary superconductor the current of the Cooper pairs satisfies the London equation. 
If the vacuum of gluodynamics behaves as a dual superconductor, then it is natural to assume 
that the abelian monopole currents should satisfy the dual London equation in the presence 
of the dual string 

E - S 2 curl j m = ® m 5(x x )5(x y )n z , (28) 

where E is the abelian electric field, 5" 1 = m P h is the dual photon mass, $ m = — is the 
electric flux of the dual string and g m is the magnetic charge of the monopole. The string 
is placed along the oz axis, the vector n z is parallel to the direction of the string and we 
assume for simplicity that the core of the string is a delta-function. 

Indeed, as shown numerically in Ref. the dual London equation is satisfied in 



the MaA projection of SU(2) gluodynamics. A recent detailed investigation |27], [| of the 
electric field profiles and the distribution of monopole currents around the string shows that 
the structure of the chromo-electric string in the MaA projection is very similar to that of 
the Abrikosov string in the superconductor. 

The following physical question is relevant: what kind of superconductor do we have. 
If niph < m mon , then the superconductor is of the second type, the Abrikosov vortices are 
attracted to each other. If m p h > m mon , then the superconductor is of the first type and 
the vortices repel each other. The computation |26| of the dual photon mass from the dual 
London equation shows that m P h ~ mmon- This means that the vacuum of gluodynamics is 
close to the border between type-I and type-II dual superconductors. The same conclusion 
is also obtained in Ref. [p8~| . A recent detailed study of the abelian flux tube |§ shows that 
the dual photon mass is definitely smaller then the mass of the monopole, and therefore the 
vacuum of SU{2) is the type-II dual superconductor (see also Ref. [|T4|| ) . 



3.3 Fact 3: Monopole Condensate 

If the vacuum of SU (2) gluodynamics in the abelian projection is similar to the dual super- 
conductor, then the value of the monopole condensate should depend on the temperature as 
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a disorder parameter: at low temperatures it should be nonzero, and it should vanish above 
the deconfmement phase transition. 

The behaviour of the monopole condensate can be studied with the help of the monopole 
creation operator. In gluodynamics this operator can be derived by means of the Frohlich 
and Marchetti construction [^] for the compact electrodynamics. For SU{2) lattice gluody- 
namics in the MA projection, it is convenient to study the effective constraint potential for 
the monopole creation operator (similar calculations were performed for compact electrody- 
namics in Ref. |J13|): 



Veff(*) = ~ H< S{* - ^ E $ mon(x)) >) 



(29) 



where ip is the monopole field. If this potential has a Higgs-like form 



V eff ($) oc A(|$| 2 - $ 



2\2 
c) i 



\,$ 2 C > o 



(30) 



then a monopole condensate exists. The dual superconductor picture predicts this behaviour 
of the effective constraint potential in the confinement phase. In the deconfmement phase 
the following form of the potential is expected: 



V eff {Q) ocm 2 |$| 2 
(no monopole condensate). 



Ai$r + ... 



m , A > ; 



(31) 
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Figure 10: V($) for the confinement (a) and the deconfinement (b) phases, Ref. [|3T| 



The numerical calculation of the effective constraint potential ( P9|) is very time-consuming, 
and therefore, in Refs. [^, |ST|] the following quantity V($) has been studied: 

(32) 



K($) = -1b(< 6($ - $ mon (x)) >). 
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Figure 11: The dependence of (a) $ c on the spatial size of the lattice for three values of (3 
and (b) of <$>f f on /3, Ref. gl[. 



Numerical calculations of this quantity were performed on lattices of size 4 • L 3 , for L = 
8, 10, 12, 14, 16 with anti-periodic boundary conditions in space directions for the abelian 
fields. Periodic boundary conditions are forbidden, since a single magnetic charge cannot 
physically exist in a closed volume with periodic boundary conditions due to the Gauss law. 
The results on finite lattices have been extrapolated to the infinite volume, since near the 
phase transition finite volume effects are very strong. 

In Figures [It] (a,b) the (right-hand side of the) effective potential ( |52"D is shown for the 
confinement and the deconfinement phases, the calculations being performed on a 4 • 12 3 
lattice. 

In the confinement phase (Figure [10] (a)), the minimum of V(<&) is shifted from zero, 
while in the deconfinement phase, the minimum is at the zero value of the monopole field 
$. The value $ c of the monopole field at which the potential has a minimum is equal to the 
value of the monopole condensate. 

The potential shown in Figure [It] (b) corresponds to a trivial potential with a minimum 
at a zero value of the field: $ c = 0. The dependence of $ c , the minimum of the potential, 
on the spatial size L of the lattice is shown in Figure |TT|(a). 

We fit the data for $ c by the formula $ c = AL a + $* n ^, where A, a and are the 
fitting parameters. It occurs that a = — 1 within statistical errors. Figure |Tl"|(b) shows the 
dependence on j3 of the value of the monopole condensate extrapolated to the infinite spatial 
volume It is clearly seen that $^ vanishes at the point of the phase transition and it 

plays the role of the order parameter. 

The monopole creation operator |32| in the monopole current representation is studied 
in Ref. [Q. First the monopole action is reconstructed from the monopole currents in the 
MA projection, and after that the expectation value of the monopole creation operator is 
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calculated in the quantum theory of monopole currents. Again, the monopole creation op- 
erator depends on the temperature as the disorder parameter. A slightly different monopole 
creation operator was studied in Refs. |34], ^5 . 



3.4 Abelian Monopole Action: Analytical examples 

In the next section we explain how to study the action of the monopoles extracted from 
SU(2) fields in the MaA projection. Now we give several examples of monopole actions 
corresponding to abelian gauge theories. 



Example 1: AD compact abelian electrodynamics 

We have already discussed the duality transformation of AD compact electrodynamics. 
There is another exact transformation which represents the partition function as a sum over 
monopole currents. This transformation was initially performed by Berezinsky |36| and by 
Kosterlitz and Thouless for the 2D XY model. Accordingly, we call it the BKT transfor- 
mation. For compact electrodynamics with the Villain action this transformation was found 
by Banks, Myerson and Kogut fl38| . The partition function for compact electrodynamics 
with the Villain action is 

Z Villam = I ' V6 Yl exp{-||d# + 27rm|| 2 }. (33) 

m(C 2 )e2Z 

As explained in Appendix C, the BKT transformation of this partition function has the form 

Z=[J VAexp{-(3\\dA\\ 2 }]x[ £ exp{-Air 2 (3 A _1 *j)}] . (34) 

-oo *j(*c 1 )e^ 

S*j=0 

Here the partition function for the non-compact gauge field A is inessential because it is 
Gaussian. All the dynamics is in the partition function for the monopole currents *j which 
are lying on the dual lattice and form closed loops (5*j = 0). 

It is possible to find the BKT transformation for compact electrodynamics with a gen- 
eral form for the action f3l| . In this general case the monopole action is non-local (see 
Appendix C). 

Example 2: Abelian Higgs theory in the London limit 



The partition function for the Abelian Higgs model is (cf., eqs. fl22"| , |23" 

+ 00 +7T 

Z AH = J VB J Vip expl-^lldEH 2 -71^ + 5 + 2™! | 2 }, (35) 

-oo -7T m(C 2 )e^ 



c We assume the Villain form for the interaction of the gauge and Higgs fields. 
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where B is the non-compact gauge field and if is the phase of the Higgs field (the radial part 
of the Higgs field is frozen, since we consider the London limit). 



After the BKT transformation fl39| , f40f , the partition function takes the form 
Z AB = const.- £ exp{-5^(j)}, S^(j) = A^j) + JL 



i(Ci)e* 

<5j=0 



(36) 



where A is the inverse lattice Laplacian and j is the current of the Higgs particles. 



Example 3: Abelian Higgs theory near the London limit 



As pointed out by T. Suzuki [141, the numerical data show that in the MaA projection 
the currents in lattice SU(2) gluodynamics behave as the currents of the Higgs particles in 
the Abelian Higgs model near the London limit. This means that the coefficient A of the 
Higgs potential A(|$| 2 — 1) , eq. (|2"3"D , is large but finite. In this case it is impossible to get an 
explicit expression for the monopole action, but there exists a 1/A expansion of this action 



Sn 



l 



J\\ +a 2 



E(E 



* -2 
J x 



X fjt=l 



where 



E 

n=l 



f, 



(n) 



A' 



(37) 



fj^ are the coefficients which can be calculated [14]. In the limit A — > +oo, this action is 
reduced to 5^f n (*j), eq.©. 



3.5 Monopole Action from "Inverse Monte— Carlo" 

The 1/A expansion of the monopole action has an important application in lattice gluody- 
namics. It is possible to show that the monopole currents in the MaA projection of SU (2) 
gluodynamics are in a sense equivalent to the currents generated by the theory with the 
action ([FJ]). This means that (at least at large distances) gluodynamics is equivalent to the 
Abelian Higgs model. The details can be found in the lecture of T. Suzuki ||14j| . Below we 
briefly describe some of these results. 

It occurs that from a given distribution of currents on the 4D lattice it is possible to 
find the action of the currents. This can be done by the Swendsen ("inverse Monte-Carlo") 
41| , [42] method. By the usual Monte-Carlo method we get an ensemble of the fields {$} 
with the probability distribution e~ s ^\ The "inverse Monte-Carlo" allows us to reconstruct 
the action S'($) from the distribution of the fields {$}. 
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The procedure of the reconstruction of the monopole action for lattice gluodynamics 
is the following. First the SU(2) gauge fields are generated by the usual Monte-Carlo 
method. Then, the MaA gauge fixing is performed and the abelian monopole currents are 
extracted from the abelian gauge fields. Finally, the inverse Monte-Carlo method is applied 



to the ensemble of these currents and the coefficients /3, 7, a 1; a 2 , a 3 , . . . of the action (|37| ) 
are determined. It occurs that at large distances the distribution of the currents is well 
described by the action (|37|). In that sense, the effective action of gluodynamics is the 
Abelian Higgs model, the monopoles play the role of the Higgs particles. In contrast to 
compact electrodynamics, the monopole potential is not infinitely deep, see Figure |l2^(a,b). 





(a) (b) 

Figure 12: Schematic representation of the monopole potential A(|$| 2 — 1) for compact 
electrodynamics (a), A — > +00; and for lattice SU(2) gluodynamics (b), A is finite. 



3.6 Various Representations for the Partition Function 

There are several equivalent representations of the partition function: 



Z J = exp{-S mm (j)} 

Sj=0 



(38) 



where S mon is given by (|3 



Representation 1 as the dual abelian Higgs model: 

As we have already said, the model (|38D is related to the (dual) Abelian Higgs theory by 
the inverse BKT transformation 



Z j oc Z AH = J VBV<5>exp{-S AH (B,<5>)} 



(39) 



where S AH is given by eq. (p3 
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Representation 2 in terms of the Nielsen-Olesen strings: 

Using the BKT transformation for the partition function ([39|), one can show that model 
|) is equivalent |43| to the string theory: 



Zi oc Z str = ]T e" 5 ™ (<T) (40) 
o-(c 2 )e^ 

<5cr=0 

S£? in9 (v) = 4 7 r 7 ( ( x ) (A + m 2 )-V)+0(i), m 2 = 1 . (41) 
Here a is the closed world sheet of the Nielsen-Olesen string. In the continuum limit [44| we 



have 

S^ mg oc J da^{x)V^{x - ~x>)da» v {x>) , (42) 
where T>$ is the free scalar propagator, (A + m 2 )!)^ (x) = 5^(x). 

Representation 3 in terms of the gauge field 9: 

Using the inverse dual transformation for the partition function (|39|) , one can get the 
representation of the monopole partition function in terms of the compact gauge field 9 
which is dual to the gauge field B, (|39|): 



Z^otZ 6 = f V9e~ s °^- (43) 

S° = — ^||d^ + 27rn|| 2 -— V(d^ + 27rn,A(d0 + 27rn)) + O(i). (44) 



Representation 4 in terms of hypergauge fields: 

Applying the duality transformation to the partition function fl43[) , we get the represen- 
tation in terms of the hypergauge (Kalb-Ramond) field 

Z j oc Z HG = Yl fvAexp{-S{A,h)} 

S(A,h) = ^||^ir + ^||^ + /.|| 2 + 0(i), (45) 

where h is the hypergauge field (h^ u (x)) interacting with the gauge field A (A^(x)); the field 
h is dual to the monopole field $, and the field A is dual to the field B. 

The action S(A, h) is invariant under gauge transformations: A — > A + da, h — * h, and 
under hypergauge transformations A ^ A — 7, h — > h + d 7 . 
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Representation 5 by fourier transformation of the string world sheets: 
This representation is intermediate between the representations fl40|) and fl39"l) : 

ZcxZ h = [VZ[V{ exp \2iri(£, a) -§[£,{) +0(\)\ , (46) 

S[tA = ^H£ + dCI| 2 + ^l|d£|| 2 , (47) 

where a are the world sheets of the Abrikosov-Nielsen-Olesen strings. Again, there exists a 
hypergauge symmetry (£ — > £ + d7, ( — > £ — 7), which reflects the fact that the string worlds 
sheets are closed (or equivalently, conservation of the magnetic flux). 



It is possible to derive a lot of physical consequences from these representations [JJj . One 
of the most interesting is that the classical string tension which is calculated from the string 
action flUD coincides, within statistical errors, with the string tension in SU (2) lattice gauge 
theory. The coefficients (3 and 7 in eq.(ffl]) are found by the inverse Monte-Carlo method, 
as we have just explained. 

3.7 Abelian Monopoles as Physical Objects 

The abelian monopoles arise in the continuum theory [l(J from singular gauge transformation 
(@)-(0) and it is not clear whether these monopoles are "real" objects. A physical object 
is something which carries action and below we only discuss the question if there are any 
correlations between abelian monopole currents and the SU(2) action. In Ref. 13] it was 



found that the total action of the SU(2) fields is correlated with the total length of the 
monopole currents, so there exists a global correlation. We now discuss the local correlations 



between the action density and the monopole currents ||5 

In lattice calculations, the monopole current j^x) lies on the dual lattice and it is natural 
to consider the correlator of the current and the dual action density^: 

1... /' . 1 .. . V ... , 1 



C x =< -Tr [j^x) ~^F a/3 (x)J > - < jfc) >< -TrF^(x) > . (48) 
The density of j^{x) strongly depends on (3 = p-, therefore it is convenient to normalize C\. 

< |Tr (j M (x) \e vlia(3 F a(5 {x)f > _ 

2 " <i^)><^W> ' 1 j 

For the static monopole we have jo(x) 7^ 0, ji(x) = 0, i = 1, 2, 3, and the fact that C 2 7^ 
means that the magnetic part S m of the SU(2) action is correlated with j^x). 

The correlator C2 is related to the relative excess of the action carried by the monopole 
current. The expectation value of S m on the monopole current is 

1/1 \ 2 

S m =< — Tr ( n M (x) -e ulxa/3 F al3 (x) \ > , (50) 



J Here and below, formulae correspond to the continuum limit implied by the lattice regularization. 
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Figure 13: The relative excess 77, eq. (|5T|) , for 10 4 lattice (a) and for 12 3 • 4 lattice (b) in the 



MA (circles) and the Polyakov abelian (squares) gauges, Ref. [45|. 



where n^(x) = j fJi (x)/\j tl (x)\. The relative excess of the action on the monopole current is 

V = ^j^, S =< ±TrF%{x) > . (51) 

where S =< ^TrF^Jx) > is the expectation value of the action, the coefficient ^ cor- 
responds to the lattice definition of the "plaquette action" < S >. If j^x) = 0, ±1, then 
T] — Ci- Since in lattice calculations at sufficiently large values of (3 the probability of 
j^,(x) = ±2 is small in the MaA projection, we have rj ~ C2 at large values of (3. From 
numerical calculations we have found that rj = C<i with 5% accuracy for j3 > 1.5 on lattices 
of sizes 10 4 and 12 3 • 4. The lattice definition of S m is 




where the summation is over the plaquettes P which are the faces of the cube C v (x); the 
cube C„(x) is dual to n v (x); Up is the plaquette matrix. Thus, S m is the average action on 
the plaquettes closest to the magnetic current (see Figure ^|(b)). The lattice definition of 
S is standard: S =< (l — ^TrUp^j >. We use these definitions of S and S m in our lattice 
calculations. 

It is interesting to study the quantity ( pT] ) in various gauges. In Figure |TB|(a) the relative 
excess of the magnetic action density near the monopole current r\ is shown for a 10 4 lattice. 
Circles correspond to the MaA projection, and squares to the Polyakov gauge. The data for 
the F\i projection coincide, within statistical errors, with those for the Polyakov gauge. 

The quantity fl5lD at finite temperatures has also been studied. In Figure |i"3"Kb) the 
relative excess of the magnetic action density near the monopole current 77 is shown for 
12 3 • 4 lattice. A gain, the circles correspond to the MaA projection and the squares to the 
Polyakov gauge. 
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Thus we have shown that in the MaA projection the abelian monopole currents are 
surrounded by regions with a high nonabelian action. This fact presumably means that the 
monopoles in the MaA projection are physical objects. It does not mean that they have to 
be real objects in the Minkowsky space. What we have found is that these currents carry 
SU(2) action in Euclidean space. It is important to understand what is the general class of 
configurations of SU(2) fields which generate the monopole currents. Some specific examples 
are known. These are instantons [46| and the BPS-monopoles (periodic instantons) [fiCfl. 



3.8 Monopoles are Dyons 

A dyon is an object which has both electric and magnetic charge. In the field of a single 
instanton the monopole currents in the MaA projection are accompanied by electric cur- 
rents [47|. The qualitative explanation of this fact is simple. Consider the (anti) self-dual 



configuration 

F^{A) = ±*F^(A) . (53) 

The MaA projection is defined Jll[] by the minimization of the off-diagonal components of the 
non-abelian gauge field, so that in the MaA gauge one can expect the abelian component of 
the commutator term l/2Tr(a 3 [A ll , A u ]) = £ A" A v to be small compared with the abelian 
field-strength f^A) = d^A^. Therefore, in the MaA projection eq. (|53D yields 



Due to eq. 



ftMu(A) « ±*f (U/ (A) . 
the monopole currents have to be correlated with electric ones, since 
f = d v f» v {A)~±dSf»M)=j™. 



(54) 



(55) 



It is known that the (anti-) instantons produce abelian monopole currents p6[ . The 
abelian monopole trajectories may go through the center of the instanton (Figure [H|(a)) or 
form a circle around it (Figure |i~4|(b)). 




W (b) 
Figure 14: The instantons (spheres) inducing monopoles (thin lines). 
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Due to eg. Q55D , the monopole currents are accompanied by electric currents in the self- 
dual fields. Therefore, the monopoles are dyons for an instanton background. However, the 
real vacuum is not an ensemble of instantons, and below we discuss the correlation of j e and 
j m in the vacuum of lattice gluodynamics. 

There is a lot of vacuum models: the instanton gas, instanton liquid, toron liquid, etc. 
Suppose that the vacuum is a "topological medium" , and there are self-dual and anti-self-dual 
regions (domains). In this vacuum the electric and the magnetic currents should correlate 
with each other. But the sign of the correlator depends on the topological charge of the 
domain. 



We define the electric current as 47 



W = ^ EM*) - - *>)]• (se) 

In the continuum limit, this definition corresponds to the usual one: = d v f^ u . The 
electric currents are conserved (d^J^ = 0) and are attached to the links of the original 
lattice. Electric currents are not quantized. 

In order to calculate the correlators of the type <J*(x)J™(x)...>, one has to define the 
electric current on the dual lattice or the magnetic current on the original lattice. We define 
the electric current on the dual lattice in the following way: 



W = Tr £ [ftx) + Jftx - ft] . (57) 



16 X £*C(y,n) 

Here the summation on the r.h.s. is over eight vertices x of the 3-dimensional cube *C(y, n) 
to which the current Jt(y) is dual. The point y lies on the dual lattice and the points x on 
the original one. 

For the topological charge density operator we use the simplest definition: 

Q(*) = W2 £ e^-^Trp^^U^ix)], (58) 

A n Ml,---,M4=-4 

where E^t 1)W is the plaquette matrix. On the dual lattice the topological charge density 
corresponding to the monopole current J™(y) is defined by taking the average over the eight 
sites nearest to the current J™(y): Q(y) = hY,Q( x )- 

X 

The simplest (connected) correlator of electric and magnetic currents is given by 

< J™J* > =< J™ J* > - < J™ >< J* > =< J™ J* >, here we used the fact that 

< J™ >=< Jf, >= due to the Lorentz invariance. The connected correlator < J™Jf, > is 
equal to zero, since J™ and have opposite parities. A scalar quantity can be constructed 
if we multiply J^Jt with the topological charge density. The corresponding connected 
correlator 

< j™j;q > c =< j™j;q > (59) 

is nonzero for the vacuum consisting of (anti-)self-dual domains (cf. eq. fl55|)). Note that we 
derived eq.fl59"D using the equalities < J" 1 >=< J e >=< Q >= 0. 
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We consider SU(2) lattice gauge theory on a 8 4 lattice with the Wilson action. We calcu- 
late HU the correlator <J™(y)J^(y)Q(y)>, using 100 statistically independent configurations 



at each value of (3. 

This correlator strongly depends on (3 and it is convenient to normalize it by dividing by 
p m p e . Here p m and p e are the monopole and the electric current densities: 



_ ' I T m ( e ) I 

V is the lattice volume (the total number of sites). 



(60) 



a-{J m J e Q)/p m p e Q-(J m J e q)/pmPe 
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Figure 15: Correlators <J m J e Q>/ p m p e and <J m J e q>/ p m p e as functions of f3, Ref. pE8 



The correlators <J m J e Q>/p m p e and <J m J e q>/p m p e (q(y) = Q(y)/\Q(y)\) are rep- 
resented in Figure [15|. As one can see from Figure |1^, the product of the electric and the 
magnetic currents is correlated with the topological charge. 

Thus our results show that in the vacuum of lattice gluodynamics the magnetic current 
is correlated with the electric current, the abelian monopoles have an electric charge. The 
sign of the electric charge depends on the sign of the topological charge density. 



4 Conclusions 



Now we briefly summarize the properties of the abelian monopole currents in the MaA 
projection of lattice SU(2) gluodynamics: 
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• Monopoles are responsible for pa 90% of the SU(2) string tension. 

• Monopole currents satisfy the London equation for a superconductor. 

• Monopoles are condensed in the confinement phase. 

• The effective monopole Lagrangian is similar to the Lagrangian of the dual Abelian- 
Higgs model. 

• Monopoles carry the SU(2) action. 

• Monopoles are dyons. 



The main conclusion which can be obtained from these facts is that the vacuum of lattice 
gluodynamics behaves as a dual superconductor: the monopole currents are condensed and 
they are responsible for confinement of color. 

We have to note that the approach discussed above is not unique and there are several 
other approaches to the confinement problem in non-abelian gauge theories. We cannot 
discuss all these here, but mention the description of confinement in terms of Z2 vortices 
49fl , the description of the QCD vacuum in terms of the non-abelian dual lagrangians [50 



and the study of QCD by means of the cumulant expansion 51 
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6 Appendix A: Differential Forms on the Lattice 

Here we briefly summarize the main notions from the theory of differential forms on the 
lattice. The calculus of differential forms was developed for field theories on the lattice in ref. 



52| . Since all lattice formulas have a direct continuum interpretation, the standard formalism 
of differential geometry can have wide applications in lattice theories. The advantage of the 
calculus of differential forms consists in the general character of the obtained expressions. 
Most of the transformations depend neither on the space-time dimension nor on the rank of 
the fields. With minor modifications, the transformations are valid for lattices of any form 
(triangular, hypercubic, random, etc). 

A differential form of rank k on the lattice is a skew-symmetrical function 0^, defined on 
k- dimensional cells of the lattice. The scalar lattice field is a 0-form. The U(l) gauge 
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field is a 1-form. The exterior differential operator d is defined as follows: 



(d0(CW!)= E ^ c 'k) (A.l) 
c fce9C fc+1 

Here <9Cfc is the boundary of the fc-cell C^. Thus, the operator d increases the rank of the 
form by one. For instance, in Wilson's U(l) theory with the dynamic variables U = e , the 
action cos(d#) depends on the plaquette d6 ) (C 2 ) constructed from the links 6{C\). 

The dual lattice is defined as follows. The sites of the dual lattice are placed at the centers 
of the D-dimensional cells of the original lattice. The object dual to an oriented fc-cell (Cfc) is 
an oriented (D— fc)-cell (*Cfc), which lies on the dual lattice and has an intersection with 
The dual of the cubic lattice is also a cubic lattice obtained by shifting the lattice along each 
axis by 1/2 of the lattice spacing. Every fc-form on the lattice corresponds to a (D— &)-form 
on the dual lattice: *4>(*Ck) = 4>{Ck)- The co-derivative is defined as follows: 

5 = *d*. (A.2) 

This operator decreases the rank of the form 

(<ty)(C»_i) = £ 4>{C k ). (A.3) 
all c k :C k _ 1 edc k 

The square of d and 5 is equal to zero: 

d 2 = 0, 

S 2 = 0. (A.4) 

The first equality is a consequence of the well-known geometrical fact: "the boundary of the 
boundary is the empty set" ; the second equality follows from the first one and the definition 

Now we discuss the lattice version of the Laplace operator, which is defined as 

A = Sd + dS. (A.5) 

This operator acts on 0-forms (scalar fields) in the same way as the usual finite-difference 
version of the continuous Laplacian. For the forms of non-zero rank the relation ( |A.5| ) is the 
generalization of the usual Laplacian. The obvious properties of 5 are 

SA = A5 

dA = Ad (A.6) 
1 = 5A _1 d + dA^ 

The last relation, called the Hodge identity, implies the widely used decomposition formula 
for an arbitrary fc-form: 

= 5 A- 1 (d0) + dA" 1 (5(f)) + fa, (A.7) 
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where <ph is a harmonic form, 5<ph = d(j)h = 0. The number of the harmonic forms depends 
on the topology of the space-time; for instance, the number of the harmonic 1-forms is equal 
to 0,1 and D respectively, for the space-time topology R D , R 0-1 <g> S 1 and S D . 

For two forms of the same rank the scalar product is introduced in a natural way: 

(0,VO= E 0(W(C fc ) (A.8) 
all c k 

This scalar product agrees with the definitions of the d and 5 operators in the sense that the 
following formula of "integration by parts" is valid: 

(d0,^) = (0,^). (A.9) 
The norm of a /c-form is defined as usual by 

M 2 = £0(cy 2 = (0,0) (a.io) 

We illustrate the above definitions by four simple examples. 

Example 1. Let us show that the general action for the compact gauge field is a periodic 
function of the angle 9p corresponding to the plaquette. By definition, 9 P = &6 + 2nk ) where 
the integer valued 2-form k is defined in such a way that —tt < Op < n. Under a gauge 
transformation we have 9' = 9 + dx and 9' p = d9 + d 2 x + ^k' = d9 + 2nk', where k' is an 
integer chosen in such a way that —7t<9'<tt. Therefore, the gauge invariance requires the 
periodicity condition on the action: S(9p) = S(9p + 2nk). 

Example 2. Note, that the definition of the norm (|A.10| ) allows us to write in concise 
form the action of the lattice theory, e.g., S = Y^c 2 [d(^4(Ci))] 2 = ||dA|| 2 for noncompact 
electrodynamics. We often use a similar notation for the Villain action. 

Example 3. Let us show that in standard compact lattice electrodynamics there exists 
a conservation law. As we have mentioned before, the compact character of the fields implies 
the existence of monopoles. The monopole charge inside an elementary 3-dimensional cube 
is defined by m = -^d9p = |d/c, where 9p = d9 + 2nk, as in Example 1. In the continuum 

limit, the above definition corresponds to the Gauss law: m = § s HdS, where S is the 
surface of the elementary cube. If j = *m, then 5j = *d**m = ^d^A; = 0. Since m is a rank 
3 form, it follows that j is a rank 1 form and the equation 

8j = (A.ll) 

means that for each site the incoming current j is equal to the outgoing one. Therefore, the 
monopole current j is conserved. It is easy to prove that j is gauge invariant. This result 



is well known for 4-dimensional lattice QED [38]. Note that all transformations considered 



above are valid for any space-time dimension D and for the field 9 of any rank k\ the current 
j in this case is a (D — k — 2)-form. For the XY model in D = 3 we have k = and the 
conservation law (|A.11|) for the 1-form j implies that the excitations, called vortices, form 
closed loops. In the 4-dimensional XY model the conserved quantity j is represented by a 
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2-form, which means that there exist excitations forming closed surfaces. These objects are 
related to "global" cosmic strings, forming closed surfaces |53j in 4-dimensional space-time. 



Example 4. To perform the BKT transformation (see Appendix C) we have to solve 
the so-called cohomological equation 

dl = n, (A. 12) 

where n is a given /c-form and I is a (k — l)-form to be found. Using the Hodge decomposition, 
we can easily show that l(n) = A~ l 5n is a particular solution which, being added to the 
general solution of the homogeneous equation, yields the general solution of ( |A.12|) : I = 
A _1 <5n + dm (m is an arbitrary (k — 2)-form). 



If n is an integer valued form, then the solution I of ( |A.12| ) can also be found in terms 



of integer numbers. We explicitly construct such a solution for k = 2. The generalization 
for arbitrary values of k is obvious. Let us assign the value to the links which belong to 
a given maximal tree on the lattice (a maximal tree is a maximal set of links which do not 
contain closed loops). To each link that does not belong to the tree we attach the value 
equal to the value of the plaquette formed by the link and the tree (there exists one and only 
one such plaquette). It is easy to see that we have thus constructed the required particular 
solution l(n). The general solution has the form / = l(n) + dm , m e For simplicity, we 
have neglected finite volume effects (harmonic forms) in the above construction. 

7 Appendix B: Duality transformation 

In this Appendix we perform the duality transformation for U(l) compact gauge theory 
with an arbitrary form of the action. Initially the duality transformation was discussed by 
Kramers and Wannier |54J] for the 2-dimensional Ising model. 

Let us consider the U(l) gauge theory with an arbitrary periodic action 5[d^], £>[..., Xp+ 
2n, . . .] — S[. . . , Xp, . . .], where 9i is the compact U(l) link field and P denotes any plaque- 
tte of the original lattice. We use the formalism of differential forms on the lattice (see 
Appendix A). We start from the partition function of the compact U(l) theory: 

Z = J V6exp{-S[d6]}. (B.l) 

— 7T 

The Fourier expansion of the function e~ s ^ yields 

Z = const. ■ JVO £ e sm e <kAe) ? (R2) 



— TT 



fc(c 2 )e^ 

where k is an integer valued two-form, and the action S d is 



S d (k) 



const. ■ J VXexp{-S[X) -i(k,X)}. (B.3) 
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The integration over the field 9 in eq.( |B.2[ ) gives the constraint 8k = 0, which can be resolved 
by introducing the new 3-form n: k = 5n. For simplicity, it is assumed here that we are 
dealing with a lattice having trivial topology (e.g. TZ A ). In the case of a lattice with a 
nontrivial topology, arbitrary harmonic forms must be added to the r.h.s. of this equation. 
Finally, changing the summation J2k(c 2 )e% — > S*n(*c 1 )e^j we S e ^ the dual representation of 
compact U(l) gauge theory: 

ZocZ d = J2 e- sd( * dn) . (B.4) 

*n(*C 1 )G^ 

Consider now compact U(l) gauge theory with the action in the Villain form eq. (|33"D . 
Now the integral is Gaussian and the dual action is S d (*dn) = \\*dn\\ 2 /4(3. 



8 Appendix C: The BKT transformation 

In this Appendix we show how to perform the BKT transformation for compact U(l) gauge 
theory with an arbitrary form for the action ( P-l|) . Let us insert the unity 



1 



J VG5{G-k) 



(here G is a real-valued two-form) into the sum ( B.2|) 

Z= ]v6exp{-S P (dO)} = const. - jvO J VG Y 5{G - k)e~ sd{G) e i{GAe) . (C.l) 



-7T — OO 



fc(c 2 )e^ 



Using the Poisson summation formula 



Y 5(G-k)= Y e2ni(G ' k) > ( c - 2 ) 



we get 

+ 7T +00 

Z = const. • J V9 J VG Y e ' Sd[G) e^ Ae+2 ^ . (C.3) 
-7r -oo fc(c 2 )e^ 

Now we perform the BKT transformation with respect to the integer valued 2-form k: 

k = m\j) + dq , dm[j] = j , d j = , (C.4) 

where q and j are one- and three-forms, respectively. First, we change the summation 
variable, Y, = X! a2 ■ Using the Hodge-de-Rahm decomposition, we adsorb 

k(c 2 )ez; i(c 3 )G^g(Ci)e^ 

dj=0 

the d-closed part of the 2-form k into the compact variable 9: 

d6 + 2vr/c = d0 n . c . + 2it5br 1 j , 6 n . c . = 9 + 2nA- l 5m[j} + 2nq . (C.5) 
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Substituting eq. flCT^ ) in eq.( p.3| ) and integrating over the noncompact field 9 n _ c ,, we get 
the following representation of the partition function: 

+00 

Z = const. • / VG e ~ Sd(G) exp{27ri(G, SA^j)} 5{5G) . (C.6) 

-00 *j{*C x )eZ 
S*j=0 

The constraint 5G = can be solved by G = 5H, where if is a real valued 3-form. Sub- 
stituting this solution into eq.( |C.6|) we obtain the final expression for the BKT-transformed 
action on the dual lattice: 

Z = const. • e' Smm( -" j) . (C.7) 



*j(*C!)e^ 

5*j=0 



where 



SmonCj) = ~^yj VHe- sd ^ ex P {2m(*H, . (C.8) 

We have used the relation d5A~ 1 j = j, Vj : dj = 0. Therefore, for the general U(l) action 

S[d6] the monopole action ( |C.8|) is nonlocal and it is expressed through the integral over the 

+00 

entire lattice ( / DH). 



As an example we consider the Villain form of the U(l) action (|33|) . Repeating all the 
steps we get the following monopole action: 

S mon Cj)=A7T 2 p(*j,A- 1 *j). (C.9) 
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